Non-singular modified gravity unifying inflation with late-time acceleration and 
universality of viscous ratio bound in F(R) theory 
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The review of basic cosmological properties of four-dimensional _F(_R)-gravity, including FRW 
equations of motion and its accelerating solutions, generalized fluid and scalar-tensor representation 
of the theory is done. Cosmological reconstruction equation is written and conditions for stability of 
cosmological solution are discussed. The overview of realistic _F(i?)-models unifying inflation with 
dark energy epoch is made. The avoidance of finite-time future singularities in such theories via 
the introduction of _R 2 -term is studied. New realistic non-singular _F(_R)-gravity unifying early-time 
inflation with late-time acceleration is presented. The exit from inflationary era in such model 
may be caused by the gravitational scenario. It is demonstrated that five-dimensional F(i?)-gravity 
considered as non-perturbative stringy effective action leads to universal relation for viscous bound 
r-\ • ratio. 
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I. INTRODUCTION 

It is widely accepted that the universe evolution passed at least two accelerating epochs: the early-time inflation 
O : and current late-time acceleration. The striking similarity between these two accelerating eras indicates that their 
origin is caused by the same fundamental scenario. The most natural and elegant explanation of the known universe 
i-C ■ evolution in unified manner may be given within modified gravity approach. It does not request any extra dark 
components like inflaton, dark energy and dark matter which should be implemented into the theory. The role of 
these dark components is played by the gravitational theory which changes its functional form at different curvature 
' scales. For instance, at the early universe the gravitational theory contains higher derivatives terms which cause 
the early-time inflation. During radiation/matter dominance higher derivatives terms are negligible if compare with 
standard General Relativity. Meanwhile, current dark energy epoch is due to the gravitational terms which dominate 
^ |* \ at very low curvature. The corresponding review on modified gravities unifying the early-time inflation and late-time 
acceleration is given in ref. [l[ . The easiest version of modified gravity which naturally unifies the early-time inflation 
and late-time acceleration is F(i?)-theory [H, 0. The number of versions of -F(i?)-gravity is known to be consistent 
with astrophysical data and local tests (for review, see Q). 

In the present paper we study different models of -F(i?)-gravity which may unify the early-time inflation with late- 
time acceleration in the consistent way. Moreover, we present new non-singular models for such unification. The 
paper is organized as follows. In the next section the general properties of _F(i?)-gravity are reviewed. Spatially-flat 
FRW equations of motion are presented. These FRW equations are rewritten in the form of standard FRW equations 
for General Relativity with generalized F(R)-Ruid. Simple accelerating cosmologies are briefly discussed. Section 
IIIII is devoted to the review of scalar-tensor representation of F(i?)-gravity (so-called Einstein frame). It is shown 
how to bound the (non-physical) anti-gravity regime in scalar-tensor representation of the theory. In section IIVI we 
propose general strategy for cosmological reconstruction of _F(i?)-gravity, i.e. derivation of its functional form which 
contains given cosmological solution. Section [V] is devoted to the study of stability of given cosmological solution. 
It gives the cosmological solution conditions which define its stability/instability. In section I VII we discuss viable 
_F(i?)-gravities unifying inflation with dark energy. The conditions for derivation of such models are given. The 
review of the known realistic unified F(i?)-theories is made. It is shown that some of these models may contain the 
finite-time future singularity. Nevertheless, the addition of i? 2 -term which supports the early-time inflation cures the 
future singularities in accord with the original proposal of refs. [^4fij|- New non-singular viable -F(i?)-gravity unifying 
the inflation with dark energy is proposed and investigated. It is shown that such theory has the instable inflationary 
era. The decay of the inflationary era is caused by the structure of gravitational theory. Section I VIII is devoted to 
the study of different, five-dimensional application of _F(i?)-gravity. Here, we discuss five-dimensional _F(i?)-gravity as 
some non-perturbative stringy gravity in frames of AdS/CFT correspondence. It is proved that for such theory like 
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for five-dimensional General Relativity the viscous ratio bound is universal. Some summary and outlook are given in 
last section. 



II. GENERAL PROPERTIES OF F(i?)-GRAVITY 

Let us briefly review general properties of four-dimensional F(R)- gravity which is known to be realistic candidate 
for the unification of early-time inflation with late-time acceleration (for review, see [if). The starting action is chosen 
to be: 



s F {R) = J d x^-g y-i^r + £« 



(i) 



For general F(i?)-gravity, one can define an effective equation of state (EoS) parameter. The spatially-flat FRW 
equations in Einstein gravity coupled with perfect fluid are: 



3 „2 1 

, p = ; 

K K 

For modified gravity, one may define an effective EoS parameter as follows: 



P=—H\ p = -—(3H 2 + 2H). (2) 



1 2H (*\ 
w eS = -1 - ^ . (3) 



The equation of motion for modified i 7 '(i?)-gravity is given by 
1 



2 9^F(R) - R^F'(R) - g^aF'(R) + V p V^(iZ) = -yT matter/1!/ . (4) 



Assuming spatially-flat FRW universe, 

i=l,2,3 

the FRW-like equations are given by 
3 „ 9 1 



ds 2 = -dt 2 + a{t) 2 ( dx% ) 2 > ( 5 ) 



1 ( 1 



p cff = ^H 2 , PcS = [3H 2 + 2H ) , (6) 



Pcff 



2 f(R) + Z(H 2 + H) f'(R) - 18 (AH 2 H + HH^j f"(R)j + p matter , (7) 
Pcff = ^ (\f( R ) ~ ( 3i/2 + h) f'(R) + 6 (sH 2 H + AH 2 + 6HH + #) f"(R) + 36 (AHH + f"'{R) 



+Pmatter ■ (8) 

There may be found several (often exact) solutions of ([6]). Nevertheless, due to presence of higher derivatives and 
non-linear terms, the structure of exact solutions is much more complicated than in General Relativity. Without any 
matter, assuming that the Ricci tensor is covariantly constant, that is, R^ v oc g^ wi Eq. Q reduces to the algebraic 
equation Q: 

= 2F(R) - RF'(R) . (9) 

If Eq. © has a solution, the Schwarzschild (or Kerr) - (anti-)de Sitter or de Sitter space is an exact vacuum solution. 
When F(R) behaves as F(R) oc R m and there is no matter, there appears the following solution: 

(m-l)(2m-l) 

H , (10) 

which gives the following effective EoS parameter: 

6m 2 — 7 m — 1 

W cf f = - T/ -TTT TT ■ (11) 

3(m — l)(2m — 1) 
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When F(R) oc R m again but if the matter with a constant EoS parameter w is included, one may get the following 
solution: 



2 in 



and the effective EoS parameter is given by 



H ^^+H i (12) 



w + 1 , 
u; eff = -l + . (13) 



This shows that modified gravity may describe early /late-time universe acceleration which could reproduce the 
quintessence-like or phantom-like behaviour. Realistic non-linear F(R) models are discussed below. 



III. SCALAR-TENSOR DESCRIPTION OF _F(_R)-GRAVITY 

It is well-known that one can rewrite -F(i?)-gravity action in the scalar-tensor form. Introducing the auxiliary field 
A, we rewrite the action (JTJ) of the F(i?)-gravity as: 

S =W-J d±x^—g{F'{A) (R-A)+ F(A)} . (14) 



2k 

By the variation over A, one obtains A = R. Substituting A = R into the action (TT41 . one can reproduce the action 
(fTj). Furthermore, we rescale the metric in the following way (conformal transformation): 

g^^e a g^, a = -lnF'(A). (15) 

Hence, the Einstein frame action is obtained: 

S e = A / d^x^—g (r - \gP°d p ad a o - V(a] 



V{a) = e'g (e-) - e^/ (g ( e -)) = _ . (16) 

Here 3 (e -<T ) is given by solving the equation a = — In (1 + /'(A)) = — InF'(A) as A = g (e _<T ). Due to the scale 
transformation (fl~5|) . there appears a coupling of the scalar field a with usual matter when matter is included. The 
mass of a is given by 

2 = 3cPF(ct) = 3 (_A AF[A)_ 1 1 

a ~2 da 2 2\F'(A) {F'(A)f F"{A) J ' 1 j 

Unless m CT is very large, there appears the large correction to the Newton law. Naively, one expects the order of the 
mass Too- could be that of the Hubble rate, that is, m a ~ H ~ 1CP 33 eV, which is very light and the correction could 
be very large. 

In ref. [1] (see also ref. @), a "realistic" F(R) model has been proposed. It may be shown however, that the model 
has a instability where the large curvature can be easily produced. This is the manifestation of the finite-time future 
singularity. In the model @, a parameter m ~ 10~ 33 cV with mass dimension is included. The parameter m plays a 
role of the effective cosmological constant. When the curvature R is large enough compared with m 2 , R ^> m 2 , F(R) 
Q has the following approximate form: 

777 2n+2 

F(R)^R-c 1 m 2 + ^ 1 ^ + 0{R- 2n ) . (18) 

Here ci, C2, and n are positive dimensionless constants. Then the potential V{a) (1 16[) has the following asymptotic 
form: 

n«>~^. (is) 

Hence, the infinite curvature R = A — > 00 corresponds to small value of the potential and therefore the large curvature 
can be easily produced. 
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Let us assume that when R is large, F(R) behaves as 

F{R) ~ F R e . (20) 

Here -Fb and e are positive constant. One also assumes e > 1 so that this term dominates if compared with General 
Relativity. Then the potential V(cr) (|16|) behaves as 

Therefore if 1 < e < 2, the potential V(a) diverges when R — > oo and therefore the large curvature is not realized so 
easily. When e = 2, V(a) takes finite value 1/Fq when R — > oo. As long as 1/Fq is large enough, the large curvature 
could be prevented. 

Note that there appears the anti-gravity regime when F'(R) is negative, which follows from Eq. (fT4"|) of ref. @ 
where first FR) model unifying inflation with dark energy was proposed. Hence, one should require 

F'(R)>0. (22) 

As long as the condition (|2"2"j) is satisfied there occurs usual gravity regime although there might appear extra force 
when rv? a (fl~7|) is small. Note also that 

Zg-gffii-AFW + iM)- (23) 

Therefore if 

= -AF'(A) + 2F(A) , (24) 

the scalar field a is on the local maximum or local minimum of the potential and therefore a can be a constant. The 
condition (f24| is nothing but the condition (j9|) for the existence of the de Sitter solution [7| . When the condition (|24|) 
is satisfied, the mass (fTT)) can be rewritten as 



2F'(A) V F"(A) 

Then when the condition (|22p for the exclusion of the anti-gravity is satisfied, the mass squared m 2 a is positive and 
therefore the scalar field is on the local minimum if 

On the other hand, if 

- A + J F$)<°- < 27 > 

the scalar field is on the local maximum of the potential and the mass squared rr? a is negative. As we will see 
later in ([551) . the condition (f2l)|) is nothing but the condition for the stability of the de Sitter space. Thus, scalar- 
tensor representation (Einstein frame) for F(i?)-gravity is constructed. Despite the mathematical equivalence of 
Jordan and Einstein frames pictures they lead to the theories which are not physically equivalent (for corresponding 
discussion, see [13]). (Or, more exactly, if physics in one frame is matched with observations and describe the 
observable accelerating universe then another frame physics becomes extremely unconventional. For instance, instead 
of acceleration the observerver sees the decceleration, or matter couples with dilaton, etc.) Note also that it is 
evident how to transform F(i?)-gravity to the effective General Relativity with generalized F(R) fluid [11]. Of 
course, the generalized fluid is gravitational fluid and it contains higher derivatives of curvature invariants. This 
fact is often ignored in the study of cosmological perturbations in i 7 '(i?)-gravity, which leads to basically erroneous 
conclusions or not well-justified approximations where fourth-order differential equations become second order ones 
similar to standard General Relativity cosmological perturbations (see, for instance, ref. [12]). For correct treatment 
the covariant, higher-derivative cosmological perturbations theory of -F(i?)-gravity was developed in ref. fl3| . However, 
it turns out that corresponding higher-derivative equations are much more complicated than the corresponding ones 
in General Relativity. 
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IV. COSMOLOGICAL RECONSTRUCTION OF MODIFIED F(i?)-GRAVITY 

Let us demonstrate that any FRW cosmology may be realized in specific F(i?)-gravity. In other words, we propose 
the general solution for inverse problem, i.e. the cosmological reconstruction of F(i?)- gravity (for the introduction 
and review, see [HI)- Let us rewrite the FRW equation ^ by using a new variable (which is often called e- folding) 
instead of the cosmological time t, N = In The variable N is related with the redshift z by 



c - N = ^ = l + z. (28) 
a 



Since ^ = H-^ and therefore = H 2 -^^ + H^-^, one can rewrite (|6]) as (for details, see [l 

- + 3 {H 2 + HH') F'{R) - 18 Uh 3 H' + H 2 (H' f + H 3 H") F"(R) + K 2 p mattor . (29) 



Here H' = dH/dN and H" = d 2 H/dN 2 . If the matter energy density /? ma ttcr is given by a sum of the fluid densities 
with constant EoS parameter Wi, one finds 

E-3(l+ii>i) -3(l+">i) -3(1+m>;W ton\ 

i i 

Let the Hubble rate is given in terms of ./V via the function g(N) as 

H = g(N) = g(-ln(l + z)) . (31) 

Then scalar curvature takes the form: R = 6g' (N)g(N) + 12g(N) 2 , which could be solved with respect to N as 
N = N(R). Using J3DJ and (J3U), one can rewrite ^ as 

= -18 (4 5 (N (i?)) 3 g' (N (R)) + g (N (Rjf g> (N (R)) 2 + g (N (R)) 3 g" (N (R))) 

( 9 (N (R)) 2 + 9' (N (R)) g (N (R)j) ^EM - EM + K 2 £ Pl0 a Q 3 ^ e - 3 ^ N ^ , (32) 



+3 



which constitutes a differential equation for F(R), where the variable is scalar curvature R. Instead of <?, if we use 
G(N) =g(N) 2 = iJ 2 , the expression (|32"T) could be a little bit simplified: 

= -9G (N (R)) (AG' (N (R)) + G" (N (R))) ^fjf ) + ( 3G (N (R)) + \g' (N (R))^ 



dR 2 V 2 / dR 

-EM. + K 2 J2 Pi0 a^ 1+Wi) e- 3 ^ +w ^ N ^ . (33) 

i 

Note that the scalar curvature is given by R = 3G'(N) + 12G(N). Hence, when we find F(R) satisfying the differential 
equation (|32|) or (|33|) . such F(R) theory admits the solution ([31]) . The reconstructed _F(ii)-gravity realizes any given 
above cosmological solution. For instance, one can present the forms of _F(i?)-theory which describe inflation or dark 
energy or even unify the early-time inflation with late-time acceleration era. The number of corresponding examples 
may be explicitly realized. However, as a rule the corresponding theory is some complicated non-analytical function 
expressed in terms of special functions. Hence, the cosmological reconstruction turns out to be approximate one. The 
exact reconstruction of ACDM universe from pure -F(i?)-gravity often leads to (almost) General Relativity [lq . 

V. STABILITY OF THE COSMOLOGICAL SOLUTION 

In the previous section, we presented a formulation to obtain the form of F(R) to reproduce the known solution 
H = H(N). Note, however, the solution H = H(N) is not always stable. In this section, the condition of the stability 
of such solution is studied. The FRW equation ^ can be rewritten by using G = H 2 and the e-folding N, as in 



= -QG(N) (4G'(JV) + G"(N)) 



d 2 F(R) 



dR 2 



R=3G'(N) + 12G(N) 



3G(N) + -G'(N) 



3 ^ r 0\ dF ( R ) 



2 v ' dR 



R=3G'(N) + 12G(N) 



F(R) 



+ ^Pmattcr (N) . (34) 

R=3G'(N) + 12G{N) 
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Let us assume the ./V-dependence of p ma tter 
the perturbation from the solution: 



is known as in (|30l) . Let a solution of (|34l) be G — Gq(N) and consider 



G(N) = G a {N) + 5G(N) . 



(35) 



Then Eq. §M$i gives 

d 2 F{R) 



= G (N) 



dR 2 
3G (N) - 



R=3G' (N) + 12G (N) 

d 2 F(R) 



5G"{N) + { 3G (N) (4G' Q (N) + G'^(N)) 



d 3 F(R) 



dR 3 



R=3G' (N) + 12G (N) 



G' (N) 



dR 2 

12G Q (N) (4G' a (N) + G'>(N)) 



R=3G' (N) + 12G (N) _ 

d 3 F(R) 



5G'(N) 



(-4Go(N) + 2G' (N) + GZ(N)) 



dR 3 
d 2 F(R) 



R=3G' (N) + 12G (N) 



dR 2 



R.=3G' (N) + 12G (N) 



I dF(R) 
3 dR 



5G. 



(36) 



fl=3G^(Af) + 12G (iV), 



Note that in this formulation, as it is assumed that the TV-dependence of /o ma tter is known, one need not consider the 
fluctuation of p ma ttcr- If the cosmological time t is used instead of N, we need to include the fluctuation of /3 ma tter- 
Thus, the stability conditions are given by 



6(4G> (N)+G'{(N)) 
36 (4G' (N) + G%{N)) 



d 3 F(R) 



dR 3 
d 3 F(R) 



R=3G' (N) + 12G (N) 



dR 3 
dF(R) 



R=3G' (N) + 12G (N) 



G (N) G (N) dR 



R=3G' (N) + 12G (N) 




fl.=3G^(AT) + 12Go(Af)y 



R=3G' {N) + 12G (N) / 



R=3G' (N) + 12G (N) / 



-6-^>0 
6G' (N) 



(37) 



12 



G (N) 



> 0. 



(38) 



In case of de Sitter space, where H and therefore Go and R — i?o = 12Go are constants, Eq. (|3T[) becomes 6 > and 
trivially satisfied and Eq. ([38} becomes 



12G 



dF(R) 



dR 



d 2 F(R) 



R=12G 



dR 2 



-Rq 



dF(R) 



R=12G 



dR 



d 2 F(R) 



R=R„ 



dR 2 



>0, 



(39) 



R=R 



which is the standard result. 

One may consider the case that a form of F (R) admits two de Sitter solutions. If one is stable but another is 
unstable, there could be a solution that matches the unstable de Sitter solution to the stable solution. If the Hubble 
rate H of the unstable solution is much larger than the Hubble rate of the stable one, the unstable solution may 
correspond to the inflation and the stable one to the late-time acceleration. Alternatively, we may directly construct 
F(i?)-gravity which admits the transition from asymptotic de Sitter universe with the large Hubble rate to another 
asymptotic de Sitter universe with the small Hubble rate. Then if the corresponding solution of F(R)-theory satisfies 
the conditions (I37|) and (|38p . such a transition surely occurs. One should bear in mind that other scenarios for exit 
from inflationary phase (for instance, due to quantum effects) may be proposed as well. 

Let us now consider the stability of de Sitter universe satisfying the condition ([9]), or equivalently (|24|) . Eq. (J9j> can 
be rewritten as 



d fF(R)\ 
dR 



Let R = Ro is a solution of (|4"U|) . Then F(R) has the following form: 

F(R) 



R 2 



= fo + f{R) (R - Ro 



(40) 



(41) 
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Here /o is a constant, which should be positive if F(R) > and n is an integer greater or equal to two:n > 2. We 
assume the function f(R) does not vanish at R = Rq, f(Ro) ^ 0. When n = 2, it follows 

- i?0+ F^)--/o + /(^o)- (42) 
Eq. (|39|) shows that the de Sitter solution is stable if 

-/o</(A )<0. (43) 

When n > 3, one gets 

Hence, we need more detailed investigation to check the stability. One now uses the expression of mr a in (|17p and 
investigates the sign of m CT in the region R ~ i?o- The expression (1251) cannot be used since it is only valid in the 
point R — Rq. Therefore 

■j.- ^-'fW u-Hr. (45) 

Eq. (|4"5j) shows that when n is an even integer, the de Sitter solution is stable if f{Ro) < but unstable if /(Rq) > 0. 
On the other hand, when n is an odd integer, the de Sitter solution is always unstable. Note, however, that when 
f(R ) < (J(Rq) > 0), we find ml > (ml < 0) if R > R but to 2 < q ( m 2 > o) if i? < i? . Therefore when 

f(Ro) < 0, R becomes smaller but when f(Ro) > 0, R becomes larger. 

Thus, it is developed the scheme which permits to check when the unification of inflation with dark energy is 
realistic one because inflarionary stage is instable. 



VI. VIABLE F(i?)-GRAVITIES UNIFYING INFLATION WITH DARK ENERGY 
A. The known realistic _F(7?)-models unifying inflation with dark energy 

In refs. [13], [3; an d [lil; viable models of F(i?)-gravity unifying the late-time acceleration and the inflation were 
proposed. To construct these models, we have required several conditions: 

1. In order to generate the inflation, one may require 

lim f(R) = -A, . (46) 

R— >oc 

Here Aj is an effective cosmological constant at the early universe and therefore it is natural to assume Aj ^> 
(lCT 33 eV) 2 . For instance, it is natural to have A^ ~ io 20 ~ 38 (eV) 2 . 

2. In order that the current cosmic acceleration could be generated, the current /(i?)-gravity is considered to be 
a small constant, that is, 

f(R ) = -2R , /'(iJo)~0. (47) 

Here R Q is the current curvature Rq ~ (lO~ 33 eV) 2 . Note that Rq > R Q due to the contribution from matter. 
In fact, if we can regard f(Ro) as an effective cosmological constant, the effective Einstein equation gives Rq = 
Rq — K 2 T m attcr- Here T mattcr is the trace of the matter energy-momentum tensor. Note that /'(Rq) need not to 
vanish exactly. Since the time scale of one-ten billion years is considered, we only require |/'(i?o)| "C (10~ 33 cV) 4 . 
Instead of the model corresponding to (|46j) . one may consider a model which satisfies 

lim f(R) = aR m , (48) 

with a positive integer m > 1 and a constant a. In order to avoid the anti-gravity f'(R) > — 1, we find a > 
and therefore f(R) should be positive at the early universe. On the other hand, Eq. (jT7| shows that f(R) is 
negative at the present universe. Therefore f(R) should cross zero in the past. 
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3. The last condition is 



lim f(R) = , 

which means that there is a flat space-time solution. 
One typical model proposed in ref. [13] and satisfying the above conditions is 

= aR 2n - m n 
J[ ' l + 7 i? n 

Here a, ft, and 7 are positive constants and n is a positive integer. Eq. (|50p gives [lj 

l/r. 

1 / 1 \ / / ; - \ 1 

i?0 = 



(49) 



(50) 



(51) 



and therefore 



f(Ro 



-2Rn = — 1 



Pi 



1 + tL 



Then it follows 



a i~ 2RqRq 2n , j8 i~ 4i?gi?Q 2 ™i?" 1 , 7 ~ 2RqRq 2rl i?™ 



(52) 



(53) 



In the model (|50|) . the correction to the Newton law could be small since the mass m a (| 1 T|) is large and is given by 
m 2 ~ iQ-i60+i09n e y2 j n so i ar S y S t em an d TO 2 ^ ]_q- 144 + 98rl eV 2 in the air on the earth [T7|. In both cases, the 
mass m a is very large if n > 2. 

As a model corresponding to (j4"6")l , we proposed the theory [l9[ 



/(#) - -ao 
One now assumes 

and 



[ tanh 



b (R - R ) 



tanh 



b R 



Ri » Rq , ai > a , bj <C 6 ■ 



When R -> or i? < i? , R/, /(-R) behaves as 

/(«) "> - 



biRi > 1 . 



a b 



aibi 



2cosh 2 (^fa) 2 cosh 2 



i?. 



and /(0) = again. When R » i?/, it follows 
/(i?) -> -2Aj = -a I 1 + tanh ' ''" 7>> " 



aj 1 + tanh 



On the other hand, when Rq <C R <C Ri , one gets 



/(#) -> -a 



1 + tanh 



boRo 



ajbjR 



2 cosh 2 (^) 



Here, we have assumed the condition (|55[) . One also finds 



f'(R) 



biRi 



-2A = -00 



a/6/ 



a/ 1 + tanh 



60 Rc 



bjRj 



1 + tanh 



2 cosh 



2 f b {R-R ) 



2 cosh 



2 f bi(R-Ri) 



(54) 
(55) 

(56) 

(57) 
(58) 
(59) 
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which has two valleys when R ~ Rq or R ~ i?/. When R = Rq, 

f'{Ro) = -a b — — y > -a/6/ - a 6 . (60) 

2 cosh 2 

On the other hand, when R = Rj, it follows 

f'(Ri) — —a/6/ T^-^ -TT- > —a/6/ — ao6 . (61) 

2cosh 2 (M«^i) N 

Due to the condition ([22]) to avoid the anti-gravity period, one obtains 

a/6/ + a 6 < 2 . (62) 
In the solar system domain, on or inside the earth, where R^$> Ro, f(R) can be approximated by 

f(R) ~ -2A cff + 2ac- b{R - Ro) . (63) 
On the other hand, since i?o "C R <C Ri, by assuming Eq. (1551) . f(R) (|54p could be also approximated by 

/(i?) - -2A + 2ae- b °( R - R ^ , (64) 

which has the same expression, after having identified Ao = A c g- and bo = b. Hence, one may check the case (|63[) only. 
The effective mass has the following form 

e b(R-R ) 

ml ~ ' (65) 

which could be very large, that is, m 2 ~ 10 1,000 eV 2 in the solar system and m 2 <~ XO 10 ' 000 ' 000 ' 000 eV 2 in the air 
surrounding the earth, and the correction to the Newton law becomes negligible. 
One may consider another model [l7j : 



( U f? \' 2k + 1 i n2fe+l 

f(R) = ^r-^ — r- ( 66 ) 

fo + h{(R-Ro) 2k+1 +R 2 o k+1 } 

It has been shown [ItJ that for k > 10 such modified gravity passes the local tests. It also unifies the early-time 
inflation with dark energy epoch. In (|B"6"j) . Ro is current curvature Rq ~ (lO~ 33 eV) . We also require 

Here Aj is the effective cosmological constant in the inflation epoch. When R ^> Aj, f(R) (JSU) behaves as 



1 /o 

A A 2 i? 2 " +1 



/(^)~-- + 72^TT- (68) 



The trace equation, which is the trace part of ((4]) looks as: 

3Df(R) =R+ 2f(R) - Rf'(R) - K 2 T matter . (69) 

Here Tetter is the trace of the matter energy- momentum tensor. For FRW metric with flat spatial part ([5]), one finds 

R ~ (to - <)- 2/(2n+3) , (70) 

which diverges at finite future time t = to- By a similar analysis, we can show that if f(R) behaves as f(R) ~ i? e for 
large i? with a constant e, a future singularity appears if e > 2 or e < 0, which is consistent with the analysis (f2"Tj) 
via the scalar-tensor form of the action. Conversely if 2 > e > 0, the singularity does not appear. Hence, adding the 
term R 2 f(R), where liniR_ i .o f(R) = ci, liniR-^oo f(R) = c 2 , to f(R) (|66|) . the future singularity (f70|) disappears. 
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Let us observe the above situation in more detail. Now we assume 

f(R) ~F + FxiT , (71) 
when R is large. Here Fq and F\ are constants where Fq may vanish but F% ^ 0. In case of i 

1 Fl = k 
fx ' 1 11 

Under the assumption (|7T) . the trace equation (69) gives 



*b = -— , = ^ , e = - (2n + 1) . (72) 



o ei n r 6 - 1 - J ^ when e < or e = 2 , . 

^i Ujtt ~ | (2 - e) .FiJF when 6 > 1 or e ^ 2 ' {<6) 

In the FRW background with flat spatial part (5), when the Hubble rate has a singularity as 

with constants ho and /3, the scalar curvature R — 6H + 12H 2 behaves as 

(F^W when 8 > 1 
R ~{ 6J S^ when8 = l • (75) 



_6J3ho_ 



when /? < 1 



In (23) or (75). /8 > 1 case corresponds to Type I (Big Rip) singularity in [20|. 1 > 8 > to Type III, > /3 > -1 to 
Type II, and 8 < -1 but 8 ^ integer to Type IV. 

The classification of finite-time future singularities used above is given in ref . [2(| : 

• Type I ("Big Rip") : For t — > t s , a — > oo, p e g — > oo and |p e ff| ~> oo. This also includes the case of /3 e ff, p e ff being 
finite at t s . 

• Type II ("sudden") : For t — > t s , a — > a s , p e s — > p s and \p e s\ —> oo 

• Type III : For t — > t s , a — > a s , p e ff — > oo and \p c s \ — > oo 

• Type IV : For t — > t s , a — > a s , p e g — > 0, \p e s\ — > and higher derivatives of H diverge. This also includes the 
case in which p e g (p e s) or both of p e g and p e s tend to some finite values, while higher derivatives of H diverge. 

Here p c g and p g arc defined by 

PcS = \h 2 , PcS = - ^ (2H + 3H 2 ) . (76) 

By substituting (75) into (73). one finds that there are two classes of consistent solutions. The first solution is 
specified by 8 = 1 and e > 1 but e ^ 2 case, which corresponds to the Big Rip (ho > and t < to) or Big Bang 
(ho < and t > to) singularity at t = to- Another one is e < 1, and 8 — ~ € l ( e — 2) (—1 < 8 < 1) case, which 
corresponds to (70) and to the II Type future singularity. In fact, we find e = — 2ra— 1 and therefore 8+1 = — 2/(2n+3). 
We should note that when e = 2, that is, /(i?) ~ R 2 , there is no any singular solution. Therefore if we add the above 
term R 2 f(R), where lim#_>o f(R) = c\, lim.R_>. 00 f(R) = c 2 , to f(R) in (66). the added term dominates and modified 
f(R) behaves as f(R) ~ -R 2 , the future singularity (|70|) disappears. We also note that if we add i?™-term with 
n — 3, 4, 5, ■ • • , the singularity becomes (in some sense) worse since this case corresponds to e = n > 1, that is Big 
Rip case. Using the potential which appears when we transform i 7 '(_R)-gravity to scalar-tensor theory [5j, it has been 
found that the future singularity may not appear in case < e < 2. 

Thus, in order to avoid the finite-tine future singularity one has to add i? 2 -term to above viable unification models: 

/(*) = aR2n - PRn + cR\ (77) 
f(R) = -oo (tanh \ + tanh O^Y) aj (tanh f ^i^M \ + tanh + cR 2 , (78) 

f(R) = ^r-^ T + cR 2 . (79) 

fo + .h{(R-Ro) 2k+1 + Rl k+1 } 
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The addition of this term has been proposed first in ref. [4J where it was shown that in this case the Big Rip singularity 
disappears. Moreover, such term which effectively cures future singularity also supports the early-time inflation. In 
other words, adding such i? 2 -term to gravitational dark energy model may lead to emergence of inflationary phase in 
the model as it was first observed in ref. Q. In case, when model already contains the inflationary era, its dynamics 
will be changed by i? 2 -term. The investigation which shows that it cures all types of future singularity has been 
done in refs. In fact, it was realized lately that some phenomenological problems [Il[ of F(i?)-dark energy 

(like consistent description of neutron stars) may be resolved in the presence of i? 2 -term. Moreover, the traditional 
phantom/quintessence (fluid/scalar) dark energy models often bring the universe to finite-time future singularity. It 
has been demonstrated in ref. (22j that the natural prescription to cure singularity in this models is again the addition 
of i? 2 -term. In other words, curing future singularity requests from specific dark energy model to be (at least, partly) 
the modified gravity! 



B. New viable F(i?)-models 

The above analysis shows that, in order to obtain a realistic and viable model, F(i?)-gravity should satisfy the 
following conditions: 

1. When R — > 0, the Einstein gravity is recovered, that is, 

, * F(R) 1 , , 

F(R) -> R that is, -> - . (80) 

R R 

This also means that there is a flat space solution as in (|4T))) . 

2. As discussed after Eq. (|4"5l) . there appears stable de Sitter solution, which corresponds to the late-time acceler- 
ation and therefore the curvature is small R ~ R^ ~ (lO -33 eV) . This requires, when R ~ R^, 

^- = foL~ fiL (R RL? n+2 + o((R- R L ) 2n+2 ) . (81) 
Here /ql and J\l are positive constants and n is a positive integer. 



3. As also discussed after Eq. (|45|) . there appears quasi-stable de Sitter solution, which corresponds to the inflation 
in the early universe and therefore the curvature is large R ~ Ri ~ (^o 16 ~ 19 GeV) . The de Sitter space should 
not be exactly stable so that the curvature decreases very slowly. This require 

F(R) foi fu (R Rif m+1 +o((R- Rj) 2m+1 ) . (82) 



R 2 

Here /oj and fu are positive constants and m is a positive integer. 
4. Following the discussion after PU]). when R — > oo, in order to avoid the curvature singularity, it is proposed 



F(i?) ^ Jo.i? 2 that is ffll^/^. (83) 



Here /oo is a positive and sufficiently small constant. Instead of (1831) . we may take 



F(R) -+ f^R 2 -* that is E^L^!±. (84) 



Here f^b is a positive constant and < e < 1. The above condition (|83|) or (|84|) prevents both of the future 
singularity and the singularity at high density matter. 

5. As in (|22[) . in order to avoid the anti-gravity, we require 

F'(R)>0, (85) 

which is rewritten as 

d (, (F(R)\\ 2 



dR^i-W" >~R- ^ 
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*- R 



FIG. 1: The qualitative behavior of versus R for a viable model. 



6. Combining the conditions (|80|) and (|85|) . one finds 



F(R) > . 



(87) 



The conditions Q] and [2] show that there must appear unstable extra de Sitter solution at R = R e (0 < R e < Rl), see 
FiglU Since the universe evolution will stop at R = Rl since the de Sitter solution R = R L is stable, the curvature 
never becomes smaller than Rl and therefore the extra de Sitter solution is not realized. The behavior of F ffi which 
satisfies the conditions dHDJ}, {ED), (JUJ), {53}, and §7} is given in FigHJ 
An example of such F(i?)-gravity is 



F(R) 
R 2 



X m (Rj; R) 



{{X m (Rj; R) - X m (Rr, Ri)) (X m (R r , R) - X m (R r , i?i)) 2n+2 



+X m (Ri;Ri) X m (Rj; Rl) 

(2m + 1) Rf n 
(R-Rif m+1 + ti 



2n+2 



.e +3 } 2,H 



2m+l 



Here n and m are integers greater or equal to unity, n, m > 1 and R\ is a parameter related with R e by 

vlu D , (2n + 2)X(R I ;Ri)X(R I ;Ri)+X(R I ;R L ) 
X (Rr, Re) ^-^ , 

It is assumed 

< R 1 < R L < i?/ . 

X (Rj; R e ) is a monotonically decreasing function of R and in the limit R — > 0, X (Rj; R e ) behaves as 

1 



X (Rj; R e ) 



R ' 



(89) 



(90) 



(91) 
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which tells that in the limit, F(R) ([88]) is (l80l) and therefore the condition [T] is satisfied. On the other hand, when 

R — > oo, 

X(R i; R e )^ i2m +2lf r ^0, (92) 

and therefore F(R) behaves as (153")) and the condition 2] is satisfied. 
When R ~ Rj, we find the behavior of (|81|). where 

i 

/OL = |-^m #l) -X"m {Rl\ Rl) 271+2 + fm + j , 

i i 2{ zVJ 

fiL = r) m T 3 (^i ^i) X m (Ri',Rl) n + foo + f (X m (Ri; — X m (Rj; Rl)) 

(2m + l) 4 < m+1 ) {Rj (R L - J? z )} 4m ( m+1 ) 



{( J R i - J R J ) 2m+1 + fi? m+1 } 



4(m+l) 



(93) 



Then the conditionals satisfied. 

On the other hand, when R ~ i?/, we also find the behavior of (|82|) , where 



v2n+2 



>.2n+2 



/oj — |(-^m (Ri)Ri) — Xm (Ri;Ri)) (X m (Ri;Ri) — X m (Ri;Rl)) 2 

+X m (Rj;Ri) X m (Rj; Rj) 2n+2 + /^ l+3 | 
■f 11 = "^2m+2 {(-^m (Ri'i Ri) — X m (Rj; Ri)) [X m (Rj; Ri) — X m (Rj; -Rl)) 2 ™ +2 

+X m (Ri;Ri)X m (R I ;R L ) n + /, 2 n+3 j j (X m (i?/; i?/) — X m (i?j; i? £ )) z 
+ (2n + 2) (X m (JJj; Jfr) - X m (i?/; i?i)) (X m (ifr; ifr) - X m (Rj; R L )f n+1 } ■ (94) 

Hence, the condition [3] is satisfied. 

We now investigate the mass of the scalar field a in order to check if the Chameleon mechanism [23] works or not. 
For this purpose, one investigates the region R\ < Rj <^ R <^ Rj. In this region, X m (Ri; R) can be approximated 
as 

X m (Rj; R) ~ — , X m (Rj; Ri) ~ — , X m (Rj; Rj) ~ — — . (95) 

rt ill ilL 

Then F(i2) can be approximated as 



F(R) f ,fn , _ i? 1 + (2n + 2)fl L / 00 2 "- 2 J^ 2 "- 2 
i? 2 ~ /o ° + i?' /n " (2n + 3)i?!i?i n+2 ~ R 2 L n+2 ' 1 J 

In the last equation, it was assumed R\ ~ i?L- If 

we assume fooR *^ /n or 

72^3»^i" +2 ; (97) 

Joe 



one gets 

ml 



4/o 



(98) 



If we consider the region inside the earth, since 1 <? ~ 6 x 10 32 eV and 1 cm ~ (2 x 10~ 5 eV) , the density is about 
p ~ lg/cm 3 ~ 5 x 10 18 eV 4 . This shows that the magnitude of the curvature is R ~ K 2 p ~ (lCP 19 eV) 2 . In the air 
on the earth, one finds p ~ l(T 6 g/cm 3 - 10 12 eV 4 , which gives Ro ~ K 2 p ~ (10 25 eV) . In the solar system, there 
could be interstellar gas. Typically, in the interstellar gas, there is one proton (or hydrogen atom) per 1cm 3 , which 
shows p ~ l(T 5 eV 4 , Ro ~ 10" 61 eV 2 . Then, the condition (JSTJ) can be easily satisfied, for example, we can chose 
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ji^ ~ MeV 2 . Then the Compton length of the scalar field becomes very small and the correction to the Newton law 
is negligible. 

There remains the condition El Since F(R) ((88]) satisfies Eqs. ([80]) and (|83]). it is clear that Eq. ((85]) or ([86]) and 
therefore the condition [S] are satisfied when R — > or R — > oo. Since F ffi is a monotonically increasing function of 
R in the region R e < R < Rl, Eq. ([85]) or (|86[) is trivially satisfied in the region. In the region i?i < Rl <C R -C -Rj , 
since -^r^ behaves as ([TO]) , we find Eq. ([55]) or ([83]) is satisfied again. Then the condition [5] seems to be satisfied in 
all the region. 

After the inflation at R = Ri, the radiation and the matter are generated generated. If the energy densities of the 
radiation and the matter dominate compared with the contribution from f(R) = R(R) — R, that is, only the first 
term in ([7]) dominates, the radiation/matter dominated universe will be realized. The densities of the radiation and 
the matter decrease rapidly compared with the contribution from /(i?)-term in the late-time and when the curvature 
arrives at R ~ Rl, the late-time acceleration occurs. 

The condition ([51]) indicates that there appears R 2 term in the very high curvature region. The i? 2 -term will 
generate the inflation besides the inflation at R = Ri if the universe started with very high curvature R^$> Rj. Since 
the inflation due to i? 2 -term is unstable, the inflation at the very early stage would stop but when the curvature 
decreases and reaches R ~ Ri, there will occur the inflationary phase again. 

Thus, we demonstrated that number of viable -F(i?)-gravity models may explain the early-time inflation with dark 
energy epoch in unified way. 



VII. UNIVERSALITY OF VISCOUS RATIO BOUND IN 5D F(i?)-GRAVITY 

In the present section, we discuss the possible applications of 5d modified _F(i?)-gravity as non-perturbative stringy 
gravity in AdS/CFT correspondence, which has provided a gravity dual framework to analyze strongly coupled gauge 



theory. Specifically, shear viscosity in the hydrodynamical limit of the gauge theories has been well studied [241 -1 
It was found that the ratio of the the coefficient of the shear viscosity r\ and the entropy density s should be larger 
than 1/4-7T, 

2>^, (99) 

S 47T 

in the Einstein gravity. It turns out this bound is not universal. The correction to the bound co ming from the higher 
order corrections, (Einstein-Gauss-Bonnet gravity or R 3 gravity,) has been calculated in refs. [28j-[30|. In this section, 
we consider _F(i?)-gravity as a prototype non-linear model which may come from non-perturbative effects of string 
theory We show the universality of the viscous bound ratio emerges even in this case. The non-universality of such 
bound is due to account of the effects in the next order of loop expansion, but our results indicate that for non-linear 
non-perturbative stringy model this may be not true. 
Let us start from _F(i?)-gravity in D dimensions: 

I = J d D xy^g~F(R) . (100) 

Later one may choose D = 5. Here we do not consider the contribution from matter. Then from ([4]), the equation of 
the motion is given by 

= \g^F(R) - R pu F'(R) + V M V„F'(ii) - g^V 2 F'{R) . (101) 
If we assume the curvatures are covariantly constant 

R = Rq , R^u = -jjg^v , (102) 
with a constant Rq, Eq. (|101j) reduces to an algebraic equation: 

= jF(R ) - R F'(R ) , (103) 



which could be solved with respect to Rq. This shows that any vacuum solution in the Einstein gravity like 
Schwarzschild solution, Kerr solution, gravitational wave, etc. is the solution of F(i?)-gravity. Especially if Rq < 0, 
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we find Schwarzschild-anti-de Sitter (S-AdS) 

-V{r)dt 2 + V- X {r)dr 2 + r 2 ^ (dx 1 ) 2 , V(r 



D-2 



ds 2 D 



i=l 



M r 2 

r D-3 72 



(104) 



Kerr-anti-de Sitter (K-AdS), and Black string-anti de Sitter solutions. Here the length scale I is related with Rq by 

D(D-l) 



i?0 — — - 



r 2 



Eq. (|104[) shows that the horizon radius ro is given by 

r = (Z 2 Af)^ 

If the radial coordinate r is redefined as 



the metric (|104j) is rewritten as 



D 



l 2 dp 2 



l 2 p 



-dt z 



Ap 2 [1 - Ml 2 p~ 

We now expand the action (|100[) from the solution g^ v — gffi (| 103[) 

a - o (0) 4- h 

Up to the third order in h^: 



p 



(105) 



(106) 



(107) 



(108) 



(109) 



S 



d D Xyf~^g 



b (Ro) - -* (R )^ V - - I — + — 2 



2" 



~F"(i? ) jv<°> "V<°> " V - (v (0) ) 2 h 



Here 



(110) 



(in) 



and v|j°^ expresses the covariant derivative given by g^u/- Then the viscosity coefficient 77 can be read from the 
coefficient in front of (d r hij) which is contained only in the 6-th term as — (F'(i? )/4) h?" (V^) 

V = F'(R Q ). (112) 

In case of F(i?)-gravity, the entropy is given by Q 

S = 4irAF'(R ) , (113) 

which gives the entropy density s as 

S 



JO) 



,s = - = inF'(Ro) ■ 



(114) 
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Combining (|114l) with (II 12[) . one gets 

S 47T 

which is not changed from the case of the Einstein gravity. This result is completely consistent with the ones in [3l[ 
and [32]. 

As an example, we consider _R 2 -gravity with a cosmological constant 



1 f o , ( D - !)P - 2 ) 
'a 



F{R:)= 2^\ R+1 T ^ + ci?2 ^ ( 116 ) 



Here the cosmological constant A is given by 



'A 

Then (fTUgj) gives 



A=- iD -T~ 2) - (^) 



^- — p — --2 C (£>-4)i 1± v i^~m — { ] 



In the signs ± of ()1 18[) . the minus sign — reduces to the Einstein gravity case in c — > limit. When the minus sign is 
chosen, in the limit we find 

D(D-l) f cD{D-l){D-A) 

Then (fTT2"j) gives 



Ro = -^-^ (1 + - \ D _' 2)ll + O (✓) > . (119) 



1 , 1 / 2cZXD - 1) , 
^=^(1 + 2^)^(1 L^j. (120) 



1 A 40c 
'a 



^2^ 1 -7T)- ^ 



Especially when D = 5, one obtains 



This result essentially agrees with the previous ones in [27], bill I22J . Note that in [32|] , the horizon radius was shifted 
from the Einstein gravity case (c = 0) by adding the correction term (c ^ 0), via account of the shift of the length 
parameter in (|117|) corresponding to the Einstein gravity to I in (|118|) . which includes the correction coming from 
c. The shift can be found by using (|118|) and (|119|) . 



which gives, in D = 5, 

(123) 



When D = 5, Eq. (|106l) shows that ro is given by ro = l^M*. Then the horizon radius ro may be shifted as 

-9 - 

r o — > (I/Ia) 2 to by the correction from ci? term and therefore the area Ao of the horizon as Ao — > {I/Ia) 2 Ao. Since 
the quantity 77 is the density on the horizon, we may regard 77 is also shifted as 

1 / 45c 
2^ { ~ 1 



'a 



v ^{l/l A )-^^— 1--^) , (124) 



which completely agrees with the previous ones in refs. [27], Hi], H3 • We should note, however, that since M in (|104[) 
is the constant of the integration the shift of the horizon radius can be always absorbed into the redefinition of M . 
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Therefore in this section, we may fix the value of the radius, which makes difference in the expression of the viscosity 
coefficient r\ but as we have seen, the ratio of 77 and the entropy density in (| 1 15[) does not depend on the radius and 
exactly agrees with the results in [13, HH HH . This is because both of the viscosity coefficient r\ and the entropy 
density s arc densities, they shift in an identical way by the shift of the horizon. For example, if we change the 
constant of the integration from Mi to M2, r\ and s shift as 

n -> {M 2 /M 1 )i 77 , s -> (Ma/Mi) 1 s . (125) 

Then in the viscous ratio the shift is canceled with each other. Hence, we proved the emergence of the universality 
viscous ratio bound for 5d F(iZ)-gravity. This may be understood also due to the fact of mathematical equivalence 
of F(i?)-gravity with scalar-tensor theory which also satisfies to universal viscous ratio bound. Nevertheless, other 
non-perturbative stringy gravity models (like F(G) gravity [33j ) should be investigated in order to understand if the 
restoration of universality of viscous ratio bound is the common property of non-perturbative models. 

VIII. DISCUSSION 

In summary, brief introduction to cosmological aspects of -F(i?)-gravity is made. Spatially-flat FRW equations 
and its typical accelerating solutions are discussed. The generalized (gravitational) fluid as well as scalar-tensor 
description of the theory under investigation is given. The emergence of the Newton law in scalar-tensor formulation 
is described. Cosmological reconstruction technique for .F(i?)-gravity is overviewed as well as general description of 
stability condition for cosmological solution. The review of known viable i ? (i?)-gravities unifying the inflation with 
dark energy epoch is given. Some of such models may lead to finite-time future singularities which may be cured by 
the addition of i? 2 -term relevant at high curvature. It is remarkable that this term not only cures future singularity 
but also contributes to realization of inflationary era. In fact, the alternative gravity dark energy model which does 
not describe inflation may also lead to inflationary phase after addition of i? 2 -term. 

New realistic non-singular i ? (i?)-gravity unifying inflation with dark energy is introduced and studied in detail. It 
is shown that inflationary epoch is unstable and decays due to gravitational action structure. Hence, modified gravity 
provides not only the universal unification scenario but also the gravitational exit from inflation. 

Section IVII addresses different, five-dimensional application of i r (i?)-gravity as non-perturbative stringy gravity 
within AdS/CFT correspondence. It is less related with previous sections which are devoted to cosmological study 
of the four-dimensional theory. It is demonstrated here that the well-known viscous ratio bound in five-dimensional 
-F(i?)-gravity is universal unlike to the Einstein-Gauss-Bonnet gravity or stringy gravity with higher-order perturbative 
gravitational corrections. 

Modified gravity unifying inflation with late-time acceleration has passed number of different checks, related with 
theoretical considerations, local tests or observational data. It is still considered as most natural candidate for 
description of the whole universe evolution in unified way. It is remarkable that just the same models may be used for 
inflation-dark energy unified description for totally different classes of theories. For instance, such unification turns 
out to be possible in recent Hofava-Lifshitz i 7 '(i?)-gravity where unified models were introduced in ref. [34[. 
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